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The exact global energy conservation laws for the full and truncated versions of the nonlinear 
electromagnetic gyrokinetic equations in general magnetic geometry are presented. In each version, 
the relation between polarization and magnetization effects in the gyrokinetic Poisson and Ampere 
' equations and the quadratic ponderomotive gyrocenter Hamiltonian is emphasized. 
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c3 ■ I. INTRODUCTION 



Nonlinear gyrokinetic Vlasov-Maxwell theory [l| presents a powerful paradigm for the theoretical and numerical 
investigations of low-frequency turbulent phenomena in strongly magnetized plasmas. The gyrokinetic Maxwell- 
I— i' Vlasov equations describe the self-consistent coupling between the gyrocenter Vlasov distribution F (in a reduced 
phase space where the gyrocenter magnetic moment is an invariant and the gyrocenter gyroangle is ignorable) and 
^p^' low-frequency electromagnetic fluctuations Ei = — V(j)i—c~^dAi/dt and Bi = V X Ai (which satisfy the gyrokinetic 
, Maxwell's equations), in the presence of a nonuniform, time- independent background magnetic field Bq. First, the 
^ ' background magnetic field Bq = Bq bo is assumed to satisfy standard space-time-scale orderings based on guiding- 
, center Hamiltonian theory 0, S, which makes gyrokinetic theory versatile enough to be modified for special plasma 
'"q_^' geometries (e.g., edge plasmas [ll, 0-01) • Next, according to the low-frequency gyrokinetic ordering [ij, it is customary 
. to use the following approximate expressions for the perturbed electric field: Ei^ ~ V^(/)i and = — bo • V(/)i — 
O I c~^dAi\\/dt. We also use the gyrokinetic fluctuation orderings [l[ c |Eij_|/i?o ~ euth, \ ~ e |Eix|, and |Bi| ei?o, 
^ . where e ^ 1 is defined in terms of the time-scale ordering for the characteristic fluctuation frequency w ~ e|r2| 
($1 = e BQ/mc is the gyrofrequency of a particle of mass m and charge e and Vth is the thermal velocity). 

^ . 

A. Energy Conservation 

(N '■ 

^ I Gyrokinetic theory provides a reduced dynamical description of nonlinear (low-frequency) plasma turbulence that 
has important analytical and numerical advantages (see Ref. 1] for further details and references). The property 
^[[^ ^ of global energy conservation of the gyrokinetic equations, in particular, is often viewed as an important test for a 
i nonlinear gyrokinetic numerical simulation code. While the total energy of the Vlasov-Maxwell equations in particle 
phase space is easily separated into the total kinetic energy of the particles and the electromagnetic field energy, 
, the explicit decomposition of the total energy of reduced Vlasov-Maxwell equations can be far less obvious. This is 
Qv^ ' because the treatment of the polarization and magnetization associated with the dynamical reduction of the Vlasov 
I equation is, in general, complicated Fortunately, the situation is not hopeless since various important sets of 
reduced Vlasov-Maxwell equations can be derived from reduced variational principles j8l-[To|. Hence, through the 
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Noether method exact conservation laws can be derived for each set of reduced Vlasov-Maxwell equations, even 
I though these equations are themselves approximate. 
5— I ' The nonlinear gyrokinetic Vlasov-Poisson [T2| - [T3 and Vlasov-Maxwell [H, [l^ equations can be numerically 
5^ , implemented by following either a "full-/" versionflTl [isj or a truncated ("5/") version fl^, H^]. The purposes of 
the present work is (i) to show that each version is associated with a numerical implementation that possesses an 
exact gyrokinetic energy conservation law and (ii) to highlight the role played by the gyrocenter polarization and 



magnetization in ensuring exact energy conservation. 



B. organization 

The remainder of the paper is organized as follows. In Sec. [Hi we present a brief review of standard Hamiltonian 
gyrokinetic theory [l| , where the gyrocenter Hamiltonian dynamics is expressed in terms of the unperturbed guiding- 
center Poisson bracket and the gyrocenter Hamiltonian (whose first and second-order terms represent the effects of 
the fluctuating fields). In Sec. IIIIl we introduce explicit expressions for the first-order gyrocenter polarization and 
magnetization associated with the gyrocenter transformation. In particular, we show how these terms appear in 
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the second-order gyrocenter Hamiltonian defined in Sec. |TT1 In Sec. IIVI we briefly review the variational derivation 
of the full-/ version of the nonlinear gyrokinetic equations and discuss the connection between the second-order 
ponderomotive term in the gyrocenter Hamiltonian and the property of exact energy conservation. First, we introduce 
the gyrocenter kinetic energy K^y = (T~y i?ogc), defined as the gyroangle-averaged part of the gyrocenter push- 
forward Tgy'^ of the (unperturbed) guiding-center kinetic energy -ffogc- Next, we present a new proof of exact energy 
conservation that emphasizes the role played by the second-order gyrocenter Hamiltonian. 

In Sec.|Vl we present the variational derivation of the truncated version of the gyrokinetic Vlasov-Maxwell equations 
(in which the second-order gyrocenter Hamiltonian is omitted from the truncated gyrokinetic Vlasov equation) , while 
the gyrokinetic Maxwell's equations retain the gyrocenter polarization and magnetization effects. From this truncated 
variational principle, we obtain an exact energy conservation law (by Noether method), which is also proved explicitly. 
Lastly, we summarize our work and discuss applications in Sec. I VII 

II. HAMILTONIAN GYROCENTER THEORY 

In this Section, we review the standard nonlinear Hamiltonian gyrokinetic theory, as derived by phase-space- 
Lagrangian Lie-transform perturbation method (additional details and references can be found in Ref. [l|). Before 
we begin with this review, we make a few remarks concerning the general theory of Hamiltonian dynamics in eight- 
dimensional extended phase space plj |. where the coordinates = {x,p;w,t) include the energy w and time t as 
canonically-conjugate coordinates. 

The extended Hamilton's equations are derived from the variational principle S J (F— "H dcr) ~ 0, where a denotes the 
Hamiltonian orbit parameter and the extended phase-space Lagrangian F — H dcr is defined in terms of the symplectic 
one- form F = F^ dz"^ (summation over repeated indices is henceforth implied) and the extended Hamiltonian Ti = 
H{x, p, t) — w. From this variational principle, we obtain the Euler-Lagrange equations uiab dz^ /da = dH/dz^, where 
the extended Lagrange (antisymmetric) matrix uJab = dTb/dz'^ — dVa/dz^ is constructed from the symplectic one- 
form F. Note that the Hamiltonian dynamics takes place on the surface H = 0, so that the energy coordinate is 
defined to be equal to the regular Hamiltonian: w = H{-k., p, t). The inverse Lagrange matrix (assuming it is regular) 
yields the extended Poisson matrix J""^ = z^}, in terms of which we express the extended Hamilton's equations 
dz"" /da = {z", T-L] — J""^ dH/dz^. Lastly, we note that an arbitrary exact exterior derivative dS may be added to 
the phase-space Lagrangian F — H dcr without changing the equations of motion. 

When perturbations are introduced in the extended phase-space Lagrangian 

F - Hda = (Fo + eFi + ••■) - (^Hq + eHi + •■•-«;) da, (1) 

both the Poisson-bracket structure { , } and the Hamiltonian H are perturbed. In the Hamiltonian formulation of 
Lie-transform phase-space Lagrangian perturbation theory, we search for a phase-space transformation z — > z (which 
induces a push- forward operation T^^) such that the new symplectic one- form F = T^^F + d5 = Fq contains no 
perturbation trace, while all perturbation effects are included in the new Hamiltonian H = T^^H ~ "Ho + eHi + 
+ • ■ • , where the second-order ponderomotive Hamiltonian 'H2 contains nonlinear polarization and magnetization 
effects [2l[ (generated by the reduced displacement T^^x). 

A. Unperturbed guiding-center Hamiltonian dynamics 

Modern gyrokinetic theory [l| is based on a sequence of two near-identity phase-space transformations. The first 
transformation involves the time-independent transformation (in eight-dimensional extended phase space) from parti- 
cle coordinates z° = {x,p;w,t) to guiding-center (gc) coordinates z^^, = (Rgc,p||gc,/^gc,^gc; W^gc*)- Here, the extended 
guiding-center Hamiltonian and the guiding-center symplectic one-form, which are respectively defined as 

T^ogc = -f^ogc - = /igc Bo +p|g^/2m - VFgc 

\ (2) 

Togc = [(e/c)Ao+p||gcbo]-dRgc + /igc(Bo/ll)d6'gc - Wgcdi J 

are independent of the guiding-center gyroangle 6gc, and thus the guiding-center magnetic moment figc (or the gy- 
roaction figcBo/D, canonically conjugate to 9gc) is an adiabatic invariant for the guiding-center motion of a charged 
particle in the time-independent magnetic field Bq = V X Aq. In addition, the particle's energy coordinate w and 
the guiding-center's energy coordinate Wgc are equal (because this transformation is time-independent). 
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The guiding-center transformation introduces the push- forward operator T~^, which transforms functions / on 
particle phase space to functions /gc = '^gc f on guiding-center phase space, such that the scalar-invariance property 
f{z) = /gc(^gc) is satisfied. Here, we introduce two important apphcations of the guiding-center push-forward '^gc- 
First, we introduce the guiding-center gyroradius 

Pgc = Tgc^x - Rgc = Pogc + £b Pigc + ■ • • , (3) 

which is expanded in powers of the guiding-center small parameter ee = Pth/^B ^ 1 (defined in terms of the 
thermal gyroradius pth = fth/l^^l and the characteristic nonuniformity length scale Lb of the background magnetic 
field). In Eq. (|3]), the lowest-order gyroradius ppg^, depends explicitly on the guiding-center gyroangle 9gc, while 
the first-order correction p^^^ has gyroangle-independent and gyroangle-dependent parts H^l- The guiding-center 
transformation, therefore, introduces an electric-dipole TTgc = e (Pgc) and an intrinsic (M) magnetic-dipole moment 
/Xgf = (e/2c) (pgj. X (Ucpgc/dt) into the guiding-center Maxwell's equations. The lowest-order electric and magnetic 
dipole moments are [23| 



TTigc \ / e(pig^) \ / ebo X (vigc)/!^ 

Pogc/ V (Pogc X 9pogc/5V> / V -Mgcbo 



(4) 



where (vigc) = Rgc± denotes the standard guiding-center magnetic-drift velocity. Note that the intrinsic magnetic- 

— M 

Ogc 

magnetic-dipole moment 



dipole moment (|4]) satisfies the variational derivation /i.^^. = — dHogc/d'Bo. In addition, the first-order guiding-center 



e 



/^igc = - ( ( Pigc X n ) + (pig,) X Rgell ) (5) 



^gc 

includes an intrinsic-magnetic-dipole contribution (involving the gyroangle-dependent part p^gc) and a moving- 

electric-dipole contribution (involving the gyroangle-independent part (Pigc)), with Rgc|| = P\\gc bo/ni. 

Next, we introduce the particle guiding-center velocity defined as the guiding-center push- forward of the particle 
velocity ^ 

Vgc = Tg/v = {Rgc + Pge, Hgclgc = Vogc + CB Vigc + • • • , (6) 

where Vogc = Rgc|| + ^f^Pogc/^^gc and Vigc — Rgc_L + vigc. We note that the particle guiding-center velocity dH) 
satisfies the relation 

Tgc^ (y Ivp) ^ y IvgcP = A^scSo + ^ + 0(4), (7) 

i.e., the guiding-center kinetic energy is identical to the particle kinetic energy (and thus the energy coordinates are 
identical: Wgc = w). 



B. Perturbed guiding-center Hamiltonian dynamics 

When the perturbed electromagnetic potentials A'( = (0i, Ai) are introduced into the guiding-center phase-space 
Lagrangian formulation ([2|), the perturbed guiding-center Hamiltonian i?igc and the perturbed guiding-center sym- 
plectic one-form Figc are defined in Eq. ([T]) as [l6j 

iJigc = e^igc = T-^(e0i) 1 

\ (8) 
Tigc = (e/c)Aigc-d(Rgc + PgJ = Tg,i [(e/c) Ai • dx] J 

Here, the electromagnetic potentials A'^ are transformed by the guiding-center push-forward operator Tg^,^ into the 
guiding-center potentials ^igc = ('/'igc, Aigc), where 

A^g, ^ Tg-l^r EE Ar(Rgc+Pgc,t). (9) 

The appearance of gyroangle dependence in the perturbed guiding-center formulation ([S]) implies that the guiding- 
center magnetic moment figc is no longer conserved, i.e., /tgc = — e (ef2/i?o) dfpigc/ dOgc + • • ■ , where the small ordering 
parameter e denotes the strength of the field fluctuations. 
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C. Gyrocenter phase-space transformation 



In order to restore the (adiabatic) invariance of the magnetic moment fj,, a second phase-space transformation 
is introduced after the guiding-center transformation. The gyrocenter transformation involves the time-dependent 
transformation from guiding-center coordinates z^^ to the gyrocenter (gy) coordinates Z"" = (R^p^^, ii,9;W,t): 



G% 



1 dG\ 

2 ^ azL 



(10) 



which is generated by the phase-space vector-components (G" , , • • • ) . Through this transformation, the fast gy- 
romotion is once again (asymptotically) decoupled from the slow gyrocenter Hamiltonian dynamics (based on the 
time-scale separation between the fast gyromotion and the characteristic fluctuation time scale) and the gyrocenter 
magnetic moment [i is now an adiabatic invariant for the gyrocenter motion of a charged particle in the perturbed 
electromagnetic fields eEi and Bq -I- eBi. 

The gyrocenter transformation (fTO|) induces a transformation on the perturbed guiding-center phase-space La- 



grangian formulation ([8|), generated by the gyrocenter push- forward operator Tgy'^, such that the new gyrocenter 
Hamiltonian dynamics is expressed in terms of the gyrocenter Hamiltonian and the gyrocenter symplectic one-form: 



r 



gy 



' KV gc 



-^gy 



W 



(11) 



:y 



where the phase-space gauge function S^y = e 5i 5*2 H (chosen to be explicitly gyroangle-dependent at all orders 

in e) generates the canonical part of the phase-space transformation ([TU]) [see Eqs. (|12p - p^ below]. By construction, 
the gyroangle-component of the gyrocenter symplectic form Fgy is (mc/e) [lAQ and the gyrocenter Hamiltonian _ffgy 
is independent of the gyrocenter gyroangle Q. 

In the present work, we use the Hamiltonian formulation of gyrokinetic theory [l[ , where the magnetic perturbation 
term Ai is removed from the perturbed guiding-center phase-space Lagrangian formulation ([8]) , so that the gyrocenter 
symplectic one- form is simply defined as Fgy = Fogc (i-e., the gyrocenter Poisson bracket is identical to the unperturbed 
guiding-center Poisson bracket, denoted as { , }gc)- The components of the phase-space generating vector fields 
(Gi, G2, • • • ) that achieve this transformation of the symplectic one- form are [l6j 



Gi — {S*!, ^gclgc + ^ -A-igc • {Rgc -|- Pgci ^gclgcj 

G2 = {5*2, 2gj.}gc — — |Rgc-|-Pgc, 5*1! X Bigc • {Rgc -I- Pgc, 2^gc}gc, 



(12) 
(13) 



where the time-invariance condition G^ = (at all orders n > 1) implies that dSn/dw 
gyrocenter Jacobian associated with the phase-space transformation pD)) is defined as 



0. We note that the 



J. 



gy 



gc 



d, 



gc Gi 



— Jo. 



igc, 



(14) 



where the unperturbed guiding-center Jacobian is Jbgc = bo • V X (Aq -|- p||gc c bo/e) and the first-order perturbed 

guiding-center Jacobian is jTigc = bp-Vx Aigc. Within the Hamiltonian formulation of gyrokinetic theory, the 
identity J7gy = Jogc in Eq. ([Tl)) is thus guaranteed by the identity Jigc = 9a(Jogc G°), where we used the canonical 
identity 9a(Jogc {Si, Zgdgc) = 0. 

The gauge functions (5*1, 6*2, • • • ) in Eqs. ([T^ - P^ are chosen to eliminate the gyroangle-dependence in the perturbed 
guiding-center Hamiltonian (at arbitrary order in the perturbation analysis). The gyrocenter Hamiltonian in Eq. (jlip 
is thus defined (up to e^) as 



e e (V-i 



(lAigcP) 



[Sl, l/'lgc} 



gc 



(15) 



where the zeroth-order term is simply the unperturbed guiding-center Hamiltonian Hogy — 1^ Bq 
first-order (covariant) guiding-center effective potential 



Igc 



= T 



gc^ 



-A,.-) . 



»lgc 



Ai 

-^Igc ^ ) 



pf:/2m, while the 



(16) 
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is expressed in terms of the guiding-center potentials ^ and the guiding-center particle velocity ([6]). Lastly, the 
first-order function 5"! in Eq. (1151) is formally defined as the solution of the first-order equation 

^^^^ EE {Si, Hogcjgc = eV'igc, (17) 

which is approximated (to lowest order in the gyrokinetic space-time-scale ordering py]) as 

Si ^ (^^) \ipi,, ^ ^ J ^i,,de. (18) 

We note that, while the gyroangle- independent potential (V'lgc) contributes to the linear (first-order) perturbed gy- 
rocenter Hamiltonian dynamics, the gyroangle-dependent potential V'lgc contributes to the (second-order) gyrocenter 
ponderomotive Hamiltonian in Eq. (jlSp . In addition, we note that the appearance of the magnetic term Ai in Eq. (jlSp 
is a reflection of the Hamiltonian formulation adopted in the present work. We now show in the next Section that the 
ponderomotive effects (associated with V'lgc) are also related to polarization and magnetization effects that appear in 
the gyrokinetic Maxwell equations. 



III. GYROCENTER POLARIZATION AND MAGNETIZATION 



The gyrocenter push-forward operator T~y^ and the pull-back operator Tgy introduced by the gyrocenter phase-space 
transformation (|10l) play a fundamental role in modern gyrokinetic theory [1] . While the gyrocenter pull-back operator 
Tgy is involved in the lowest-order integration of the gyrokinetic Vlasov equation , the push- forward operator Tgy^ 
is intimately connected with the reduced polarization and magnetization effects [3, [2lj in the gyrokinetic Maxwell 
equations. 



A. Gyrocenter Gyroradius 



Reduced polarization and magnetization effects in nonlinear gyrokinetic theory are associated with the first-order 
gyrocenter gyroradius Pigy, defined as the gyrocenter push- forward of the particle guiding-center position 



e p 



igy 



~ gy V g 



gc 



fl_ /dpo dSi _ dpo dSi 



r 9Si cbo 

bo t; — + —fT X VSi 

op\\ eHo 



(19) 



which is expressed exclusively in terms of Si (and therefore i/'igc) and we henceforth use the lowest-order guiding-center 
gyroradius pg = ppg^, and guiding-center particle velocity Vq = Vogc- 

The first-order gyrocenter gyroradius (|19p can be decomposed into a gyroangle-independent part 



, . n d / dSi\ e d I ~ 



(20) 



and a gyroangle-dependent part 



Plgy - Pig 



(Pigy) 



dSi cbo 



eBo 



X V^i 



B^ 



d_ 



Po 



dSi 



d_ 



Po 



dSi 
W 



Po 



dSi 
W 



(21) 



We note that the gyroangle-independent part (PU)) has only perpendicular components while the gyroangle-dependent 
part (|2ip has both parallel and perpendicular components. 



B. Gyrocenter Dipole Moments 

We now make a few comments concerning the role played by the first-order gyrocenter gyroradius (|19p in modern 
gyrokinetic theory. First, according to the general reduced Vlasov-Maxwell theory [2l| . the gyrocenter transformation 
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introduces additional electric and magnetic dipole moments (for each particle species), which are defined as 

TTigy = e (pigy) = - ^gJ^{pQ V'lgc) , (22) 

Migy = e ^pigy X — ^ = - ^pigy X r! — ^ + e {p,^^) X — - 

^ f ("» - - 7 ( t ^> - I ■ 

Here, we note that the magnetic-dipole moment ()23p is naturally decomposed [251126} into an intrinsic part (involving 
Pigy) and a moving-electric-dipole part (involving (pi^y)). The gyrocenter definitions ([^ - (^5]) are therefore natural 
extensions of the guiding-center electric and magnetic dipole moments (|4]) . 

In order to get a clearer picture of these gyrocenter dipole moments, we apply the zero-Larmor-radius (ZLR) limit 
(Po 0) on the gyrocenter moments ((22|) -([23 | . In this limit, the first-order gyrocenter electric-dipole and intrinsic 
magnetic-dipole moments are TTigy ~ e p^^ and /xj^y ~ — /iBi/i?o, where the first-order gyrocenter gyroradius (in 
the ZLR limit) is 

Hence, in the ZLR gyrofluid approximation the gyrofluid polarization charge density is expressed (at lowest 
order) as 



= - V 



mn 



Eli + ^ X Bij 



en^, (25) 



where n and u\\ denote the gyrofluid density and parallel gyrofluid velocity. In Eq. (|25p. we easily recognize the 
perpendicular perturbed electric and magnetic-flutter contributions as well as the compressional contribution from 
the parallel perturbed magnetic field. On the other hand, the ZLR gyrocenter polarization velocity 

dpi^ _ c ( aEi_L , P||bo ^ dBn \ / 1 dAi±_\ ^ cbo 



dt Bq^ \ dt mc dt j \c dt J Bp ' 

includes the standard polarization velocity (with its magnetic-fiutter counter-part) and the inductive part of the Ex B 
velocity, i.e., (Ei + V0i) X cbo/Bo = ~ dAi/dt X bo/Bo. 

C. Gyrocenter Push-forward of the Guiding-center Particle Velocity 

Next, we consider the gyrocenter push-forward of the guiding-center velocity 

e Aigc ~ 



Tgy vo = vo - e f js*!, voj 



gc mc 



= {R + Po. (^ogc + ee(^ige) + •••)}^^ + e^Pigy + (26) 
where we used the relation 

with the total first-order polarization-displacement velocity [2l| defined as 

dgcPigy _ ( ^,1 dpi 



dt 



{pigy, ^ogcj ^ - ~it^ + {^igy ^^*sc}^^- (27) 
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The gyroangle average of the gyrocenter push-forward relation ([26)) is expressed as 

(Tgyivo) = {r, (ffogc + ee(^igc) + ' ' ■)} + Ky) + •••)' (28) 

which displays the characteristic separation [2l|, [2^ of the gyrocenter velocity R = Rgc + e Rigy + • • • , where 

Rlgy ^ (r, 6(^^1^0)1 = - ^^^^ bo + ^ X V(V'ige), (29) 

I. J gc m c -Do 

and the total first-order gyrocenter polarization velocity dgc (Pgy )/ dt. Note that, while the gyroangle-averaged effective 

potential (ipigc) contributes to the slow gyrocenter dynamics, the gyroangle-dependent potential V'lgc contributes to 
the gyrocenter polarization velocity (through (Pigy))- 

D. Polarization and Magnetization Effects in the Gyrocenter Ponderomotive Hamiltonian 

Lastly, we point out that polarization and magnetization effects associated with the first-order gyrocenter gyroradius 
()f 9p can also be made to appear explicitly in the gyrocenter Hamiltonian (|f 5p by rewriting the second-order gyrocenter 
Hamiltonian as 

H2sy = - ^(^Igc- {R+Po, e (V'lgc)}^j - I (plgy- (Elgc + ^ XBige)), (30) 

where we have defined the guiding-center perturbed electromagnetic fields 

Elgc = -Vc^-igc - c-^dAigc/dt ] 

\ ■ (31) 

Blgc EEE V X Algc J 

We note that, within the Hamiltonian formulation of gyrokinetic theory, the ponderomotive term in Eq. (j30p appears 
in the form of a nonlinear finite-Larmor-radius correction p9| as follows. 

First, we introduce the gyrocenter push- forward Tgy^^/^igc of the effective first-order potential (|T6| . defined as 

■^igy = TgyVigc = i/'igc - - Algc - {R + Pq, e (-iAigc)} 

C J gc 

- e Pigy • (Elgc + ^ X Blgc) - Algc • Pigy) + • • • , (32) 

where we used the definitions pip for the guiding-center electromagnetic fields. Next, by replacing the guiding- 
center perturbation potential ■0igc with the gyrocenter perturbation potential p2p in the gyrocenter Hamiltonian and 
ignoring the exact time derivative in Eq. p2p . since the Hamiltonian is defined up to an exact time derivative within 
Hamilton's variational principle, the gyrocenter Hamiltonian is now expressed as 

2 

Hgy = Hogc + ee (V-igy) + e (pig^ • (Eigc + ^ X Bigc)) , (33) 
where the second term on the right side involves the gyroangle average of the gyrocenter potential (|32p . 

IV. FULL HAMILTONIAN GYROKINETIC FORMULATION 

In this Section, we review the exact energy conservation law for the full-/ version of the nonlinear gyrokinetic 
Vlasov-Maxwell equations within the context of the gyrocenter polarization and magnetization effects introduced by 
the gyrocenter phase-space transformation (jlOl) . The explicit proof of energy conservation is presented here with a 
focus on the role played by gyrocenter polarization and magnetization effects as they appear in both the gyrocenter 
Hamiltonian and the gyrokinetic Maxwell equations. 
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A. Variational formulation 



The variational derivation of the full-/ version of the gyrokinetic Vlasov-Maxwell equations is based the gyrokinetic 
action functional [1, |^ 

Ay ^ - f -^gy ^gy d'Z + [e^ |Eip - |Bo + 6 B,\') , (34) 



where the extended gyrocenter Hamiltonian Hgy = Hgy — is defined in Eq. (ITSl) . while the extended gyrocenter 
Vlasov distribution J'gy = F 5{W — Hgy) guarantees that the energy constraint Hgy = is satisfied. 

In the full-/ version of the Hamiltonian formulation of gyrokinetic theory the time evolution of F is determined 
by the nonlinear gyrokinetic Vlasov equation 

f + ^gy}^^ = 0, (35) 
where the gyrocenter Hamiltonian (jl5|) is now written as 

-f^gy = ^ + f^Bo + ee(?Aigc) - ye (i^gyi/'igc) (36) 

where the first-order gyrocenter Lie-derivative operator, generated by the first-order field (fT2|) . is defined in terms of 
an arbitrary function g as 

£gyg={Si,g} + -Aige- {r + Po, 5} ■ (37) 

^ J gc C ^ gc 

Next, we introduce the variational expression 

SHgy = ee((5i/'igc) - e^e (i;gy(5Vigc) = ee (T^^f^i/'igc) , (38) 
where the (truncated) gyrocenter push- forward operator T^"'^ is defined as 

Tr^5 = .9 - e £gyg- (39) 
and we used the identity associated with the definition (jlSp for Si: 

'{Si, eSMgc) = {{Si, {SSi, %gc}gc}^J ^ (^{sSi, {Si, Hogc}gc}^J + ^({^1, ^Sij^c) 

[{SSi, e^Pig,}^^) , (40) 

which is obtained by using the Jacobi identity for the guiding-center Poisson bracket and omitting the exact time 
derivative from the variation i5Hgy (which is allowed in Hamiltonian dynamics). Equation ((38|) can also be used to 
write the functional variation of the gyrocenter Hamiltonian ([36]) 

SHgy = 

where the guiding-center delta function 



(T-'Sl^) SMr) - {t7' {Sl ^)) •<5Ai(r) ] d'r, (41) 



S'g, EE T-^5^(x-r) = <5-^(R + Po-r) (42) 



is evaluated at an arbitrary field point r and the identity 

5Vgc = J d^r(sMr^,t) - ^•<5Ai(r,t)) ^ 

follows from the definition (fTBl). 
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By substituting the variation (|4T|) in the variational principle derived from the gyrokinetic functional ([34)) . we obtain 
the gyrokinetic Maxwell's equations [s^l, which describe the self-consistent coupling between the gyrocenter Vlasov 
distribution F and the electromagnetic perturbation fields Ei and Bi: 

V X ,B. + , BO - 1 . - 4, f = 4. . f (T.- SI,) ) , ,44, 

where we use the gyrocenter push-forward definitions 

Tr'^^gc = ^gc + ePigy-V^^e, (45) 

and 



Tr' (vo si) ^ (vo + e Rigy) ^ -t- e (^S^ -fp,^^ + vo p,^y ■ VS^j . (46) 

When the polarization and magnetization terms (associated with Pigy) are separated from the gyrocenter source terms 
[2l| . the gyrokinetic Maxwell's equations (|43l) - (|4l|) can also be expressed as 



e V-Di = 47r / eF (S^) , (47) 
Jz 

V X (Bo + 6Hi) - ^ ^ ^ 47r 1^ e F ((^ S^,) + e il,,y (<5|,)) , (48) 

where the reduced electromagnetic fields Di = Ei + 47rPigy and Hi = Bi — 47rMigy are defined in terms of the 
first-order gyrocenter polarization and magnetization (using Sgc — S^) 

( ^igy \ 

F d^p, (49) 

\ /^igy / 

which are expressed in terms of gyrocenter phase-space moments of the first-order gyrocenter electric-dipole and 
magnetic-dipole moments (|22|) - (|231) . We note that the gyrocenter source terms in Eqs. (|47)) -(|48 | now only include the 
effects due to guiding-center polarization ((^gc)) and magnetization ((vq S^)), in addition to the perturbed gyrocenter 
current density associated with the first-order gyrocenter velocity Rigy defined in Eq. (1^^ . 




B. Gyrokinetic Energy Conservation 

The exact global gyrokinetic energy conservation law dEgy/dt = for the full- / gy rokinetic Vlasov-Maxwell equa- 
tions (|35l) and (|43l) -(lt4 l) is expressed in terms of the total gyrokinetic energy (TgI. Isif 

E,y = J^i^' + |Bo + eBip) + F (i/gy - ee(T7Vigc)), (50) 

which includes higher-order terms associated with the electric-field energy. This expression was derived by Noether 
method [8, 9] and, hence, one can explicitly show that dEgy/dt = at all orders in e. The explicit proof presented 
heer emphasizes the role played by the second-order ponderomotive gyrocenter Hamiltonian, on the one hand, and 
the gyrocenter polarization and magnetization effects in the gyrokinetic Maxwell's equations, on the other hand. 

1. Gyrocenter kinetic energy 

The Vlasov term in Eq. ([50|) introduces the gyrocenter "kinetic" energy 

^gy = (Tgy^ogc) = Hgy - ee(T7Vigc) 

Vo \ , 2 _ / ^ f , 1 



= Hogc - ee (Aigc-— ^ + e e (£gy ( 0igc - -^tpigc ) ), (51) 
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which removes the first-order contribution from 0igc and alters the second-order ponderomotive Hamiltonian. By 
rearranging the second-order terms in Eq. (j5ip . we obtain the convenient expression 



1 

2m 



po - e - Aigc 



2n 



\ ^Igc, 01gc \ - \ Slgc, 6gc \ ) , 
I. J gc 1^ J gc / 



(52) 



where we use the notation ($igc, '^igc) = J ('/'igc ^igc) d9, with ^igc — Aigc -vo/c, and we introduced the gyrocenter 
"canonical" kinetic energy 



1 

2m 



Po - e - Aigc 



PI 

2m 



ee ( Ai 



vo 

c 



2 mc^ 



(lAigcH, 



(53) 



defined in terms of the guiding-center momentum pg = mvp, which satisfies the guiding-center identity ([7]). 

The last terms in Eq. ([5^ involve separate ponderomotive electric and magnetic contributions whose physical 
significance is explained as follows. First, we return to the expression (j5ip for the gyrocenter kinetic energy and, using 
the low-frequency definition (£gy(/)igc) — (Pigc * Eigc), we obtain another form of the gyrocenter kinetic energy ([51]): 



Kgy - Hogc 



YU^Pigy 



■ Eig 



ePigy X — ) -Bi 



igc 



(54) 



Here, we simply note that the negative sign between the polarization and magnetization energies is needed in order 
to recover the standard expression for the second-order macroscopic electromagnetic energy [25| 



(Ei-Di 



87r 



|Bi| 



igy 



•El - Mw-Bi) 



(55) 



where the first-order gyrocenter polarization and magnetization are defined in Eq. (|49p . 



2. Proof of exact energy conservation 

We now proceed with the explicit proof of energy conservation dEgy/dt = 0. In contrast to the previous proof of 
energy conservation [3l| (which ignored the higher-order terms associated with Ei -I- c V0i ^ 0), the proof presented 
below focuses on the role played by quadratic nonlinearities in the gyrocenter Hamiltonian (|36|) . here expressed as 

Hgy = i/ogc + ee ^(V'lgc) - I (i^gyV'igc)^ , (56) 

where we introduce the nonlinearity ordering parameter e~in Eq. (j56p . which is either t — e (for a nonlinear gyrocenter 
Hamiltonian) or e~ = (for a truncated gyrocenter Hamiltonian) . 

By defining the push-forward operator Tr^ = 1 ^'e£gy associated with the nonlinearity parameter e~, the expression 
for dEgy/dt is separated into four energy-exchange terms: 

^ = + rH + V^ + Pa, (57) 



dt 



where the Vlasov term is 



the Hamiltonian term is 



the Poisson term is 



dF 

Vf ^ I -g^Hgy, (58) 



z 



^ I ^ - " (I (T.-Vigc) ) ) , (59) 



z 



11 



and the Ampere term is 



Va 



d-'x — 



e 9Ai 



47r dt 



V X B 



c ^ 



■ V^^s^ c 



(61) 



Note that we obtain these energy-exchange terms directly from Eq. (jSOl) by rearranging terms and integrating by 
parts (where surface terms are expected to vanish based on suitable boundary conditions). The separations shown 
in Eqs. ()58|) -()61 |) combine terms of similar nature and facihtates the proof and its interpretation. The proof that 
dEgy/dt = in Eq. (|57p uses the operator-commutation identity 



' ' dt 



- T-i ( ^9 



dt 



d 



dt 



dSi 



dt 



e dAio 



c dt 



^ J s 



(62) 



between the partial time derivative operator and the push- forward operator (acting on an arbitrary function g), which 
involves the explicit time dependence of the gyrocenter transformation. 

First, using the gyrocenter Vlasov equation ([55)1 . the Vlasov term ([55]) can be written as 



Vf = - 



[ {(FiJgy), EE 0, 



(63) 



when we use general property /^{/ifflgc = of the guiding-center Poisson bracket (for arbitrary functions / and 
g). We note that Eq. (j63p vanishes under all versions of the nonlinear gyrokinetic equations (i.e., for arbitrary e). 
Second, we express the partial time derivative of the gyrocenter Hamiltonian (|36p in the Hamiltonian term (j59p . using 
Eq. as 



dt 



e e 



f g(V'lgc) 

I dt 



e e 



^ \ dt 



and the Hamiltonian term ()59p thus becomes 

'dipigc 



Vh 



e F 



dt 



dt 



(TfVigc: 



e F 



' ' dt 



(64) 



where the operator-commutation definition (j62p is used. Here, we note that this energy-exchange term is non- 
vanishing only when the gyrocenter Hamiltonian p6p retains the second-order ponderomotive contribution it ^ 0), 
i.e., [T(^i, d/dt] = 0. Third, we use the partial time derivative of the gyrokinetic Poisson's equation (|43p to obtain 



47r dt 

so that the Poisson term (l60l) becomes 



dF 
'dt 



e F 



' ' dt 



(TrVigc) 



dF 
'dt 



F 



' ' dt 



'Igc 



((T^Vigc) - (TrVigc)) 



(65) 



Here, the first term on the right side of Eq. (1651) involves the gyrocenter polarization in the gyrokinetic Poisson's 
equation (j43p . while the remaining terms vanish for the nonlinear case e = e. Fourth, we use the gyrokinetic Ampere's 
equation (|44|) to obtain 



d^x dAi 
1^ "df 



V X B 



e 9Ei 
c "dT 



eFiJ 



\ dt ' c 



and we use the operator-commutation formula (|62p . with g = Cigc to obtain the Ampere term 



Va ^ e eF 



9^ 



(66) 



Here, the first term on the right side of Eq. (|66p involves the gyrocenter magnetization in the gyrokinetic Ampere's 
equation (j44p . while the remaining terms vanish for the nonlinear case e~= e. 
By combining Eqs. (I551) - (|66l) into Eq. ((57)) . we obtain 



dt 



eit-e) 



dF 
'dt 



gy 'rigc 



) + F 



d 



{dt 



("^gy ^igc 



F 



£ 



d_ 
'di 



V'lgc 



(67) 
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where we expanded the push-forward operator -'^ = 1 — e £gy (associated with the gyrokinetic polarization and 
magnetization) and the push- forward operator Tr^ = 1 — e £^y (associated with the nonhnear gyrocenter Hamil- 
tonian) . Equation (|57)) clearly shows that only the case of a nonlinear gyrocenter Hamiltonian (e~ = e) satisfies the 
exact gyrokinetic energy conservation law dEgy/dt = 0. The gyrokinetic energy conservation law also implies that 
the energy-transfer processes can be represented in terms of the pair-wise couplings 

^ Vh ^ Va, (68) 

where electric and magnetic fluctuations separately contribute to the evolution of the gyrocenter Hamiltonian (and 
thus to the evolution of the gyrocenter Vlasov distribution). 

Lastly, we note that if the nonlinearity parameter e~ is set to zero in the gyrocenter Hamiltonian (|56p. we obtain 
dEgy/dt = ©(e^), i.e., the non-conservation of gyrokinetic energy arises at the same order as the nonlinear time scale 
associated with the turbulent gyrokinetic dynamics, which can lead to unphysical saturation levels [1]. In particular, 
the Hamiltonian energy transfer term (j64p vanishes for t ^ 0, which interferes with the nonlinear energy-transfer 
processes ([68|l . 



V. TRUNCATED HAMILTONIAN GYROKINETIC FORMULATION 



Most numerical simulations of the nonlinear gyrokinetic Vlasov-Poisson and Vlasov-Maxwell equations ignore the 
quadratic gyrocenter Hamiltonian in the gyrokinetic Vlasov equation [i.e., = in Eq. (|56p ]. Under these conditions, 
it is important to ask whether there still exists an exact gyrokinetic energy conservation law. The answer is yes but 
only if we use an explicit truncation scheme in the variational principle (j34p for the nonlinear gyrokinetic equations. 



A. Truncated gyrokinetic Vlasov-Maxwell equations 



In the truncated (Sf) version, the gyrocenter Vlasov distribution is expressed as F = Fq + e Fi, where _Fo is the 
reference distribution function and -Fi denotes the departure from Fq driven by linear and nonlinear effects. Here, 
-Fb is assumed to be time independent and satisfies the equilibrium condition {-Fb, ^o}gc = 0, i.e., it is an arbitrary 
function, not necessarily Maxwellian, of the constants of the motion for guiding-center Hamiltonian dynamics. 

The truncated gyrokinetic (trgy) Vlasov equation is written as 







dgcFi 
dt 



{( 



Fn 



eFi 



(V'lgc)} , 
J gc 



(69) 



where dgc/dt = d/dt+ { , -ffogcjgc denotes the (unperturbed) guiding-center Vlasov operator, the linear drive is 
e {Fq, e (V'igc)}gc and the nonlinear drive 



{Fi, e(i/'igc)}g 



^xV(^rgc)-VF, -Kebo-fvF,^-^ 

Bo " V 9p\i 



V(V'lgc) , 



(70) 



which combines the perpendicular nonlinear drive and the so-called parallel nonlinearity. The perturbed gyrocenter 
Vlasov distribution Fi is coupled to the truncated gyrokinetic Maxwell equations 



eV^0i = -47r 
V X (Bo -l-eBi) = 47r 



z 
e 

z c 



(Fo-feFi) (Si) - 6Fo {£gySl) 



(Fo + eFi) (voSl) - eFo (£gy (vq <5|,) ) 



(71) 

(72) 



where the effects of gyrocenter polarization and magnetization are represented by the gyrocenter operator £gy Here, 
the terms (<5gc) and (voSg^) are associated with guiding- center polarization and magnetization, respectively, while 
(£gy(5gj.) and ("Cgy(vo5|c)) associated with gyrocenter po\a.Tization and magnetization. 

The truncated gyrokinetic Vlasov-Maxwell equations (p5|) - (l7^ are obtained from a truncated variational principle 
based on the action functional 



A 



trgy 



"87 



(e^ |V0i 



iBn 



eBil 



eFo {£gy ■01gc) 



/ -^trgy (iJogc + ee(V^igc) - d^Z, 



(73) 
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where the last term represents the truncated gyrocenter Vlasov part, with J^trgy = {Fq + ei^i) S{W — -fftrgy) denoting 
the truncated gyrocenter Vlasov distribution in extended phase space (with i?tigy = -f^ogc + e e (i/'igc)). We note 
that the term T H2gy in the full gyrokinetic action functional (|34|) is replaced with Fq -ff2gy in Eq- (ESI), and we have 
inserted the low-frequency approximation Ei = — V(/)i. The second integral term (involving only the background 
Vlasov distribution Fq), therefore, does not contribute to the variational derivation of the truncated gyrokinetic Vlasov 
equation ()69|) . Instead, it only contributes to the gyrocenter polarization and magnetization effects in the truncated 
gyrokinetic Maxwell's equations ([7T|) - ([7^ . 



B. Exact energy conservation law for truncated gyrokinetic equations 



The exact gyrokinetic energy conservation law dEtrgy/dt = for the truncated gyrokinetic Vlasov-Maxwell equations 
(p5)) - (|72p . which is derived by Noether method from the variational principle ([75| . The truncated total gyrokinetic 



energy is 



E. 



'trgy 



(Fo + eFO Lfftrgy - ee(0ige> + Fo 



' d^x 
"8^ 



IBr 



eBil 



(74) 



where the second-order gyrocenter kinetic energy is 

K2sy = i?2gy + 6 (£gy(?!)lgc) 

(lAigcP) + 



2 



2 mc^ 



2^1 



The expression for dEngy/dt is again separated into four energy-exchange terms: 

dE^^j-gy 



where 



dt 

VtrF - e 



Vt,F + Vt, 



■H 



Vtrd 



Ogc 



eFo 



ee(V'igc)) 



dFi 



trA, 



0, 



Vtr 



■H 



e Fq (^£gy 

dK2gy 



01gc j , 

t^Aigc vo 
dt ' c 



- v< 



tiA- 



(75) 

(76) 

(77) 
(78) 
(79) 
(80) 



Note that Eq. ([75]) is obtained from Eq. (|65l) by keeping only the first term while Eq. (1751) is obtained from Eq. ([TO)) 
by setting in the last term and partially cancelling the remaining terms. In Eq. (j80p . we used the relation 



2gy 



fdH; 



2gy 



dt 



V dt 

-e ( £ 



£ 



Igc 



gy " 



dt 



+ e 



( dipigc 
dt 



dt 



dSi 
'dT 

£ 



Plgc 



d_ 

di 



gc/ 



and Eqs. ([78|) - (l79|) for Vtrij) a-nd Pha- Hence, we obtain dEtigy/dt = in Eq. ([76]) and it is thus possible to construct a 
truncated version of the nonlinear gyrokinetic Vlasov-Maxwell equations that conserves energy exactly. Note that the 
energy-transfer processes ([55]) are now replaced by the truncated nonlinear processes 'Ptr<p ^ VtrH ^ 'PtrA expressed 
in terms of the background distribution Fq. 



C. Limiting forms 

As a last topic, we now discuss useful limiting forms of the truncated total gyrokinetic energy (1741) . In particular, 
we discuss the use of a truncated form of the guiding-center Poisson bracket, which leads to simplified expressions for 
the gyrocenter kinetic energy (1521) . 
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In all expressions involving the guiding-center Poisson bracket { , }gc appearing in the polarization and magneti- 
zation contributions in the gyrokinetic Maxwell equations (|7ip - (|72p . it is common to use the approximation 



({/, 3}gc) 



B 



where / and g are arbitrary functions, and the Poisson-bracket terms associated with the slower parallel and drift 
dynamics are omitted. Under this approximation, the second-order gyrocenter kinetic energy (j75p becomes 



2gy 



(81) 



where we used the identity (xigc) 
therefore becomes 



(X?gc> 



(Xigc)^ for Xigc 



trgy 



igc and ^igc. The total gyrokinetic energy ([74 
(lAigcH 



^2 ^'Fq 

2 mc^ 



2T 



^?gc 



"87 



IBn 



■F 



2T 



) - (01gc)') 
- ((4c) - (6gc) 



(82) 



where we have assumed that the background gyrocenter Vlasov distribution Fp is Maxwellian in (with temperature 
T). Similar limiting forms have been discussed elsewhere by Hahm et al. [1, [l3|- Here, we easily recognize three sets 
of terms, each set with an immediate physical interpretation. The first group of terms corresponds to the truncated 
form of the canonical kinetic energy ()53p . the second group (involving ^i) includes the electric-field energy and its 
gyrocenter polarization correction, while the third group (involving Ai) includes the magnetic- field energy and its 
gyrocenter magnetization correction. 



VI. SUMMARY AND FUTURE WORK 



In conclusion, we have presented explicit proofs of the exact energy conservation laws for the full-/ and truncated-^/ 
versions of nonlinear gyrokinetic Vlasov-Maxwell equations. These proofs relied on the intimate connection between 
the gyrocenter ponderomotive Hamiltonian and the gyrocenter polarization and magnetization effects. In Sec. IIII[ we 
introduced the first-order gyrocenter gyroradius (jl9p in terms of which the gyrocenter polarization and magnetization 
effects are introduced explicitly in the gyrocenter Hamiltonian (155]) and the gyrokinetic Maxwell equations (H7|) -(H5 | . 

In Sees. In] and IIIIl we showed that each dynamical reduction associated with unperturbed (guiding-center) and 
perturbed (gyrocenter) gyromotion introduces polarization and magnetization effects into the reduced Hamiltonian 
and the reduced Maxwell's equations. This connection between dynamical reduction and polarization-magnetization 
effects has been studied in the general context (2l| . in the high-frequency oscillation-center reduction ^J, and in the 
case of the bounce-motion reduction [s^l. In Sees. HVl and fVl we elucidated the connection between exact energy 
conservation and the accurate treatment of gyrocenter polarization and magnetization. 

Future work will look at the prospect of constructing a variational principle from which the energy-conserving 
generalization of the nonlinear Frieman-Chen equations [35| might be derived. 
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